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Glueballs and the superfluid phase of Two-Color QCD 
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Abstract. We present the first results on scalar glueballs in cold, dense matter using lattice simulations of 
two color QCD. The simulations are carried out on a 6 3 x 12 lattice and use a standard hybrid molecular 
dynamics algorithm for staggered fermions for two values of quark mass. The glueball correlators are evalu- 
ated via a multi-step smearing procedure. The amplitude of the glueball correlator peaks in correspondence 
with the zero temperature chiral transition, fi c = m^/2, and the propagators change in a significant way 
in the superfluid phase, while the Polyakov loop is mearly insensitive to the transition. Standard analysis 
suggest that lowest mass in the ++ gluonic channel decreases in the superfluid phase, but these obser- 
vations need to be confirmed on larger and more elongated lattices These results indicate that a nonzero 
density induces nontrivial modifications of the gluonic medium. 
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1 Introduction 

The study of the phases of QCD has applications to a wide 
array of subjects, ranging from heavy ion collision from 
FAIR to RHIC to LHC energies, the physics of compact 
stars, QCD phase transitions in the early universe, and 
more theoretical scenarios in which strongly coupled gauge 
theories might serve as a laboratory for physics beyond the 
standard model pQ. 

Lattice field theory provides a solid computational frame- 
work for these studies, in all those cases in which the de- 
terminant is positive and importance sampling can be ap- 
plied. Unfortunately, at high baryon density the weight of 
the functional integral is complex, which forbids the usage 
of simple importance sampling in the computation. 

At high temperature strong fluctuations of the baryon 
number makes it possible to study the physics of nonzero 
baryon density, and in this way we can explore the por- 
tion of the phase diagram which is analytically connected 
with the zero density quark gluon plasma region. Unfortu- 
nately, this excludes many different other phases of QCD, 
which appear at low temperatures and higher densities, 
see e.g. [5] for an overview. 

Two color QCD is an important exception: the color 
SU(2) group in the fundamental representation is pseudo- 
real and in this case the determinant is real even when 
/i / and standard hybrid Monte Carlo simulations are 
possible. So lattice calculations can confront and extend 
fermionic model studies [5]. 



However, very little happened in the field after the pio- 
neering studies of the 80's OIMS]- The subject was revived 
after the work of the late 90's [THE]: these authors noted 
that phenomena associated by diquark condensations at 
high density might well be much more important than 
previously though, if one goes beyond one gluon exchange 
when calculating the relevant four fermion couplings. Two 
color QCD was then reconsidered as a laboratory for a toy 
study of diquark condensation mechanisms 9 . Symme- 
tries and spectrum were analysed in a quantitative way, 
and lattice studies of diquark condensation were carried 
out rrolfTT] . Further studies include RMT analysis [12l[l"3] 
and studies of the Dirac spectrum [T3] . Moreover, the fact 
that jit c (T = 0) = makes the model amenable to a 

xPT analysis [15] . In turn, chiral perturbation theory mo- 
tivated studies produced results for the EoS and Gluon 
Condensate [IB]. The model (with adjoint fermions) has 
been also used as a laboratory for the study of new univer- 
sality classes [T7] . The phase diagram has been computed 
in the chiral limit by use a new cluster algorithm |18j . 
Model studies of the vector sector [TH] prompted studies 
of vector spectroscopy on the lattice [2Tfll2"Tir2"2"ll23| . 

In short summary, model analysis, chiral perturbation 
theory, and the feasibility of lattice calculations conspired 
to make the fermionic sector of two-color QCD reason- 
ably well understood, with numerical results in quantita- 
tive agreement with analytic calculations. 

In comparison with the fermionic sector, the gluody- 
namics of two-color QCD is much less well known. Stud- 
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ies have been carried over at high temperature, where 
the finite /i transition is similar to the one at finite T 
[25,26 , while at low temperature the behaviour of the 
Polyakov loop and of the gluon propagator indicates that 
the hadronic/superfhhd transition and the deconfining tran- 
sition are well separated phenomena [27] . 

In this study we present a quantitative analysis of the 
scalar plaquette-plaquette correlator, with the main goal 
of making progress on our understanding of the phase di- 
agram of QC2D at low temperature, and to get further 
insight into the nature of the gauge field dynamics at high 
baryon density, in a region which is distinctly different 
from the ordinary, high temperature, high baryon density 
phase. Glueballs have been considered as probes of the 
thermal medium in [28 , 29J . 

Previous work on the subject has discussed at length 
the similarities and differences between the phase dia- 
grams of two and three colors QCD. In short, while the 
properties of the matter above T c are substantially inde- 
pendent on the number of colors the nature of the low tem- 
perature, high density phase seems to depend strongly on 
the number of colors. At the same time, though, many of 
the differences appear in the fermionic sector, where differ- 
ent couplings produce different pattern of the symmetries 
realisation. The gluonic sector is not directly related with 
these different symmetries, and might well exibit charac- 
teristics independent on the number of colors. This is one 
of the motivation for a lattice study. On a more general 
grounds, two color QCD at high density affords an exam- 
ple, as we will see, of a chiral transition clearly separated 
from the deconfining one. 

This paper is organized as follows: in the next Section 
we review a few basic facts about two color QCD and its 
phase diagram; Section [3] is devoted to the details of the 
simulations, while Section 2] deals with observables and 
measurement strategies. Smearing |33j is introduced and 
reviewed there. Section [5] contains the results: we discuss 
the behaviour of the scalar gluonic correlators constructed 
by the smearing technique in order to excite glueballs and 
study their behavior as a function of baryon density. Re- 
sults for the Polyakov loop are presented as well. In Sec- 
tion [6] we summarise and discuss our findings, and outline 
future research. 

A preliminary account of some of the results have ap- 
peared in [301I3T] . 



2 Action, symmetries, phase diagram 



Let us briefly remind here a few features of the thermo- 
dynamics of the two color model on which we will rely in 
the following. The above mentioned reviews and literature 
can of course be consulted for more details. 



The lattice version of the Action reads [9]: 

Skin= E ^ftWWxMx + P)- 

x{x)U${x - v)x{x - i>)]+ 

X 

X(x)e-»u}(x-t) x (x-t)} 

where a chemical potential per quark was introduced in 
the usual way[H[3"!3]. 

At \i = the Action enjoys the lattice equivalent of 
Pauli-Giirsey symmetry: 

r 2 u fl T 2 = u; (1) 

where t 2 is Pauli matrix . 
The Action can be recast 

S kin = £ f ^[X e {x)U v (x)X {x + u)- 

x even, v— 1.3 

X e (x)Ul(x - v)X {x - P)} + 

^[^W(o e \)U t (x)X {x + t)- 

x even 

X e {x)f*^)Ut{x-t)X {x-t)] 
with the definitions 

X e = (Xe,-xtr 2 ) : X = (*° tr ). (2) 

' Mo 

which make even more transparent the equivalence of quarks 
and antiquarks in the model. 

At fj, — The (Pauli-Giirsey) symmetry implies the de- 
generacy of mesons and baryons (diquarks) with opposite 
parities: 

m n = m qq (3) 

In turn, this produces only one discernible condensate, 
< ipip > 2 + < iprp > 2 , whose direction is chosen by the 
mass term. 

Consider now the thermodynamics of the model. 

At zero density, high temperature there is a transition 
separating the ordinary hadronic phase from a quark gluon 
plasma. The characteristics of the high temperature and 
low temperature phases seem to be fairly independent on 
the number of colors, with T c of the order of a typical 
mass scale (of course the precise nature of the transition 
is not). This remains true for a nonzero chemical potential, 
till the temperature if "high enough" , most likely till the 
transition is between a hadronic phase and a quark gluon 
plasma phase. 

At zero temperature standard arguments indicate that 
the there is no Fermi sea until the quark chemical potential 
exceeds the mass of the lowest state which carries baryon 
number, i.e half the mass of the baryon which is degener- 
ate with the pion: // c = m 7r /2. For /i > [i c the attractions 
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Phase diagram of QC 2 D 
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Fig. 1. Sketchy view of the phase diagram of two color QCD: 
we recognize (at least) three phases. A Hadronic phase, a BEC 
phase, and a Quark Gluon Plasma phase. The phases are char- 
acterised by the values of the < ijjip > and < ipip > conden- 
sates, the Polyakov loop and the number density. (In this plot 
we do not distinguish between transition and crossover lines) 

in the diquark channel leads to the formation of a diquark 
condensate < ipip >. For two color QCD < ipip > is color 
neutral, like the chiral condensate itself: when p, exceeds 
p c , < iptp > rotates to < ipip >, and numerical simu- 
lations nicely follow the prediction of chiral perturbation 
theory for the chiral condensate, mass spectrum and ther- 
modynamics. 
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Fig. 2. Effective masses for the p and the it mesons in the 
normal phase (p = 0.0) and bare quark mass m, = 0.05, from 
the two different procedures labeled as el and e2, see text. 



Values fj, = 0.25,0.30,0.32,0.35,0.90 were simulated 
only at m q = 0.07, with a slightly reduced statistics of 
30000 steps. For the analysis we used a bin size of 7000, 
for which the errors level off. 

To set the scale, we have measured the fermionic spec- 
trum at p = 0.0. 

We have constructed the zero momentum propagators 
C(t) for the pion and the p mesons in the standard way, 
assumed the form 



3 The simulation 

The simulations were performed using a standard hybrid 
molecular dynamics algorithm for staggered fermions us- 
ing the same code as of ref. [SJ, which corresponds to 
Nf = 8 continuum flavors. 

Details on relevant parameters are summarised in ta- 
ble [TJ Smaller dt were needed in the supercritical region 
and trajectory lengths ranged in most cases from half to 
one. 



C{t) = Acosh[m(t- N t /2)\ 



(4) 



which is valid for euclidean times N t larger than the re- 
ciprocal of the excited masses, and used two different def- 
inition of the (time dependent) effective masses m e //(i), 
one exploiting C(t)/C(t + 1) 

C(t)/C(t l 1) - - ^/ 2 )] f5 ) 

c(tyc{t + 1) - coshK//( ^ + 1 _ Nt/2)] (5) 

the other C(t)/C(t + 2), according to 



Table 1. Simulation parameters 



lattice size 


6 3 x 12 


P 


1.3 


configurations 


40000 


sweeps 


10 


dt 


0.05 - 0.02 


thermalization 


1000 configurations 


quark mass m q 


.05, .07 


chemical potential p 


.0, .2, .25, .30, .32, .35, .4, .6, .9 


smearing weights ui 


.025, .05, .07, .1, .2, .3, .4 


smearing steps 


1,2,3,4 



C{t)/C(t + 2) 



cosh[m e/ /(f)(t- N t /2)] 
cosh[m e //(i)(i + 2 - N t /2)} 



(G) 



We show the results for the p and the pion at m q = 0.05, in 
Fig. [21 where we label the effective masses extracted by use 
of the two definitions as el and e2 respectively. The two 
procedures give comparable results on a moderate range 
of times, indicating that the oscillatory contribution for 
the opposite parity partner, if any, is very small. 

^From the plateaus in the plots we read off = 
0.56(2) and m p = 1.4(1). A similar analysis at m q = 0.07 
gives m p — 1-5(2) and = 0.64. In both cases m p and 
m T are well resolved, meaning that we are still within the 
range of chiral perturbation theory. 
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4 Glueball measurements: operators and 
smearing 

The operators commonly used for measuring scalar glu- 
onic correlators exciting glueball masses are Wilson loops. 
For simplicity we restricted ourselves to plaquette-like op- 
erators that can be built from four links. 

Simple glueball wave functions such as the plaquette 
have the disadvantage that they only have small overlaps 
with the lowest-lying glueball states, and that these over- 
laps become rapidly smaller as the lattice spacing is de- 
creased. 

Furthermore the plaquette couples strongly to ultra- 
violet fluctuations, increasing the noise in the correlators; 
to have reliable glueball correlation functions at different 
distances, it is mandatory to reduce these statistical fluc- 
tuations. 

Different smoothing procedures aimed at removing the 
unphysical short-distance fluctuations have been introduced 
a short review can be found in [37]. One of this is the 
smearing method, originally proposed for SU(3) in [33] . 
It consists in the construction of correlation functions for 
operators which are a functional of the field smeared in 
space and not in time. 

For each link U^x) of a configuration the product of 
the other three links defining a plaquette is considered, 
then these products are summed over the four choices of 
plaquettes orthogonal to the time axis; the resulting ma- 
trix is added with a weight w to the original link U^x). 
This linear combination, projected back on the gauge group, 
is assumed as new link variable (the so called "fat link"). 
The procedure can be iterated, say N s times. 

This operation both enlarges the "size" of the ordi- 
nary plaquette operator and reduces the noise associated 
with the additional degree of freedom due to the gauge 
invariance. Only spatial links participate on the averag- 
ing. Thus the transfer matrix for the smeared operators 
is unaffected and it remains positive definite, if this was 
the case. The values of the smearing coefficient w and the 
number of iterations N s are tuned in order to optimise the 
performance of the method [3"4ll35[ l36] . 

The glueball operators are defined by means of the 
plaquettes Pij(x,t) on ij plane as in the following: 



These operators are the lattice counterparts of the con- 
tinuum 0(3) <g> Z(2), J PC = 0++ and 2++ respectively. 
Hence, we can label the corresponding excitations as ++ 
and 2 ++ states. Of course this correspondence is not one 
to one but infinite to one. What is usually measured is 
the lowest excitation in the corresponding representation 
of the cubic group, corresponding to the fundamental glue- 
ball state. 

Unfortunately, but not unexpectedly, the 2 ++ results 
turned out to be too noisy to provide any significant infor- 
mation. Because of this, we will base the following analysis 
on the ++ results alone. 

With dynamical quarks, as in our case, mixing is pos- 
sible between purely gluonic and fermionic operators. So 
strictly speaking we should not talk about "glueball masses" , 
but, rather, of the lowest excitation in that particular 
chanell, be it fermionic, gluonic or a mix of the two. This 
said, for the sake of brevity, in the following we will some- 
times speak of "glueballs" , meaning what has just been 
explained. 

Following the methods previously depicted, a smeared 
operator is obtained replacing the ordinary link variables 
U{x) in the plaquette operator Pij(x,t) by the fat link 
U'(x). Higher levels of smearing are obtained by varying 
the weight w and by iterating the procedure N s times. 

Glueball masses are calculated from the Euclidean time 
behaviour of their correlation functions: 



o++ W=tr£ 



P 12 (x,t) + P 23 {x,t) + P 13 (x,t) 



(7) 



which transforms according to the (one-dimensional) 
irreducible representations of the of the relevant cubic 
point group (on a lattice the full rotational symmetry is 
broken down to only cubic symmetry), and 



e + (*)=tr]r 

X 



P 12 (x,t) - Pi 3 (x,t) 



P 12 {x,t) + P 23 {x,t) - 2P 13 {x,t) 



(8) 



(9) 



which transform both according to the E ++ (2-dimensional) 
representation. Note that all operators are projected onto 
zero momentum by averaging over all spatial sites. 



N t -1 



C 



0++ w^E(<^ 0+ >+^ 0+ >)> 



T = 



(10) 



<^° ++ (r) > 2 ) 
where we exploited the cluster property 



lim o++ (t + ;)0 o++ (t) =0° (t) 2 

t—>oo 



(11) 



as well as the translational invariance of <j>° (t). In the 
following whenever no confusion could arise we will refer 
to the propagator eqfTOl simply as C(t). 

We have analysed the glueball correlations for differ- 
ent smearing parameters. One observation which allows a 
more compact analysis, as well as offers more control on 
the results, was proposed in ref. [35]. There, it has been 
shown that the two-dimensional parameters space of the 
number of sweeps N s and the smearing weight w may 
be reduced to a single dimension via the parameter T s = 
N s x w, in a certain range < w < w max , < N s < N™ ax , 
where w max and N™ ax depend, of course, on the specific 
application under consideration. 

In our case, we have observed that the unidimcnsional 
parametrisation remains true till w < 0.3, irrespective of 
the number of smearing steps, at least within the allowed 
range of smearing steps N s < 5/2, where S is the spatial 
size of the lattice. Fig. [3] demonstrates this property in the 
normal (/i = 0.0) and in the superfluid phase (/i = 0.6) 
for the amplitude C(0) of the correlators built for our 
smallest quark mass m q — 0.05. We have observed that 
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Fig. 3. Amplitude of the ++ correlator for smearing weight w 
and smearing steps N s as a function of the smearing parameter 
T s = N s w, in the normal (fi = 0.0) and superfluid phase (/j = 
0.6), for m q = 0.05. Note the universal behaviour as a function 
of T a in the two phases. Note also that smearing affects the 
results at [i = 0.6 more than those at /i = 0.0 



this behaviour holds true also for other observables, and 
mass m q = 0.07. 

To be on the safe side, we will base our discussions on 
results within the range of validity of the unidimensional 
parametrisation: in conclusion, w < 0.3, N s < 3. 

The same Fig. [3] allows a comparison of the effects of 
smearing in the two phases: we note that the smeared 
results drop much faster as function of T a after the tran- 
sition (n — 0.6) than in the normal phase (/x = 0.0). We 
will come back to this point in the next section. 



5 Results 

We devote this Section to the presentation of the results 
for the glueballs and the Polyakov loop. 

The raw results for the plaquette correlators 

, m-i 

p w ^E< ^ 0++ ( t + ^ 0++ w > ( 12 ) 

for our smaller mass m q = 0.05 are presented in Fig. 2J 
The results in the two phases, the normal phase at /z = 0.0 
and the superfluid phase, fi = 0.4,0.6, show the expected 
long distance plateaux. 

This suggests to base our analysis on the propagators 
defined in eq. 1101 However we will also consider the pla- 
quette correlators eq. 1121 which will require leaving the 
vacuum contribution as a fitting parameter. 

The latter - fits to the plaquette correlators - will be 
our only analysis in the case of the critical propagators, 
since they show a peculiar behaviour, which we will de- 
scribe a bit more fully at the end of this Section. 

To get a first feeling about the results in the two phases 
we plot in Fig. [5] the ++ glueball correlators (eq fT0|) . for 
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Fig. 5. Propagators of the scalar glueball at mass m q = 
0.05(0.07), upper (lower) in the normal and superfluid phase, 
for smearing parameters w — 0.3, T a — 1, normalised to one 
at zero distance. Either diagrams demonstrate the decrease of 
the mass of the lowest excitation in the scalar glueball channel 
when going from the normal to the superfluid phase. 



the two masses we have inspected, smearing parameters 
as indicated, and for the two phases. Other smeared data 
behave similarly. We normalise the propagators to one at 
zero distance (the behaviour of the amplitudes will be dis- 
cussed below) to put in evidence the decay rate as a func- 
tion of euclidean time. It is immediate to notice that the 
glueball correlator decays much slower in the superfluid 
phase, indicating that the lowest excitation in the gluonic 
channel becomes lighter. 

In the remaining part of this section we discuss first 
the behaviour of the amplitudes of the correlators, then we 
analyze their long Euclidean time behaviour. The former 
indicate the critical point, the latter allow a quantitative 
estimate of the spectrum in the two phases. 



5.1 Amplitude 

Fig.[5]shows the amplitudes C(0) of the correlators, eq.fTUl 
as a function of the chemical potential, for the raw data, 



and for a representative value of T s . Other smeared data 
behave similarly. Remember that our estimates of the pion 
mass, as well as previous results [9], indicate \x c (m q = 
0.05) = 0.28(1) and fi c {m q = 0.07) = 0.32(2). We note 
that the amplitudes peak in correspondence with these 
values. 

This is consistent with their definition. In fact, C(0) is 
defined as 



1 N t -1 



> 



<0° (r)> 2 ) (13) 



hence, it is related with the plaquette susceptibility, which 
is a standard indicator of criticality. 

A second observation concerns the amplitude in the 
superfluid phase, where the amplitude itself gets smaller. 

Moreover, see again Fig. [51 which shows that the effect 
of the smearing becomes more pronounced in the super- 
fluid phase. Smearing cancels ultraviolet fluctuations and 
preserves long distance physics. Naively, any observable 
which has only short distance components would disap- 
pear after smearing. Since in the high density phase the 
residual amplitudes is comparably smaller than in the low 
density phase, one might conclude that the relative weight 
of long distance, low mass components in the high density 
phase is comparatively smaller than in the low density 
phase - in other words, if the amplitudes would become 
zero, the glueball would decouple completely. 

At even larger fi value the results increase and ap- 
proach the quenched ones, in agreement with the conclu- 
sions in ref. [25] . 

The main conclusion from the consideration of the am- 
plitudes is that the either the transition and the superfluid 
phase are clearly seen also in a pure gluonic channel. This 
demonstrates the effect of a finite density of baryons on the 
gauge fields: since our operators do not contain explicitly 
the chemical potential, this dependence is only possible 
if the gauge field themselves are affected by a density of 
baryons, both in a pure phase and at the critical point. 

These results indicate a sizable modification in the glu- 
onic sector induced by a finite density of baryons. It is then 
very natural to ask ourselves if the confining properties of 
the theory are affected as well. To address this question, 
even at a very preliminary level, we study the behaviour of 
the Polyakov loop, which is an exact order parameter for 
confinement in the quenched theory, and an indicator of 
enhanced screening properties with dynamical fermions. 

We show in Fig. [TJthe behaviour of the Polyakov Loop 
at mass m q = 0.05 as a function of fi. The effect, if any, 
is very modest and more visible for the spatial loop. As 
a comparison, we show also the spatial Polyakov loop. To 
attempt to put in evidence bubble or ordered structures 
we have also computed, and plotted in the same figure, 
the average of the absolute modulus of the local loops. 
They stay constant as well. These results suggest that the 
confining properties of the superfluid phase appearing at 
fi q > TOjt/2 are the same as those of the hadronic phase, in 
agreement with |27j . In the same paper [27] the Polyakov 
loop is shown to increase from near-zero values at even 
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Fig. 6. Amplitudes of the correlator as a function of fi for raw 
and smeared data, and quark mass m q = 0.05(0.07), upper 
(lower) diagram. The amplitudes, related with the plaquette 
susceptibilities, peak at the critical point /i c = m,r/2 and de- 
crease in the superfluid phase. These features demonstrate that 
the phase transition is visible in the gluonic sector, and that a 
high density of baryons affect the gluon mass spectrum in the 
superfluid phase. 
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Fig. 7. Polyakov loop P, and spacial Polyakov loop P^ a t mass 
m q — 0.05 as a function of /i. It is also shown the (rescaled) 
average of the local absolute values of the same observables. 
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Fig. 8. Glueball effective masses as a function of Euclidean 
time, for /i = 0.0 and m, = 0.05. Results for T s — N 3 x w = 
0,0.1,0.2,0.3,1.2 are plotted slightly displaced one from the 
other, from left to right. 



larger /i, implying that the confining property of the vac- 
uum is lost at very high density. In this paper we only con- 
sider the BEC superfluid phase appearing at [i q > 
and we do not discuss any further transition or crossover 
at larger /i. 

5.2 Glueballs in the normal phase 

The measurements of the glueballs in normal conditions 
are notoriously difficult, since the glueballs themselves are 
rather heavy. To be more specific, the plaquette correla- 
tors are basically constant within errors. To put in ev- 
idence the exponential decay one needs to subtract the 
long distance plateau. This is accomplished by exploiting 
the cluster property, and subtracting the square of the 
plaquette. Once this is done a signal for the mass - the 
value of the exponent - can be detected, while the errors 



of the fits on the amplitudes remain large. This procedure 
is justifiable only when all of the excitations are massive, 
which we expect to be the case in the pure phases, while 
in the critical region the presence of soft, massless excita- 
tions might call for some additional care. In any case, the 
validity of the procedure can be assessed a posteriori from 
the quality of the exponential fits. 

We present some results here for the sake of compari- 
son with the analysis in the superfluid phase, which is the 
main goal of this study. 

We have started with an effective mass analysis, where 
a time-dependent mass m e / / (i) is extracted according to 
eq. 

We have studied the smearing dependence of the ef- 
fective masses extracted at different Eulidean times. The 
effective mass should approach the true lowest excitation 
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le-03 



le-04 r 



le-05 



le-06 



le-07 



H = 0.0 m q = 0.05 .— b- 



Fit [1:4] 
Fit [2:4] 
Fit [3:4] 



Table 2. Results for the ++ Glueball at n = 



1 2 3 4 5 6 

Euclidean time 

Fig. 9. Details of the fits for different time interval for /i = 0.0, 
m q — 0.05, smearing weight w = .2 and N 3 — 2 (T a = 0.2). 
Fits to propagators with other smearing parameters behave 
similarly, see text for details, and Table [2] for a summary of 
the results. 



either at large T s and/or large Euclidean time, and con- 
sidering the simultaneous dependence on these parameters 
should in principle allow a good control on this measure. 

We show the results for the effective masses extracted 
from correlators at m, = 0.05, for T s = 0.0,0.1,0.2,0.3,1.2 
in Fig. [5] The results are very noisy, anyway it is clear that 
in the normal phase the glueball is heavy as expected, and 
our results allow us to place at least a lower bound for its 
value. 

Next, we performed direct fits of the data to the form 
of eq. [5] In attempt at improving the quality of the results 
we considered, along side with T s = 0.6 and 1.2, two higher 
values of T s . 

For m q = 0.05 the fits for t > 2 were satisfactory, while 
those including t = 1 were not adequate to describe the 
largest time distances, as apparent also from a higher % 2 . 
As an example, in Fig. [9]we show some fits for the param- 
eters w=0.2 and N s — 2 (T s — 0.4). The resulting masses 
where 1.32(2) for the fit performed in the time interval 
[1:4] , which is clearly missing the large time separation 
behaviour, while the fits in the interval [2:4], and [3:4] 
gives 0.68 and 0.63, in good agreement with each other, 
but errors on the fit can be hardly estimated, given the 
large errors on the propagators themselves. The behaviour 
demonstrated in Fig. [9] is generic for the m q — 0.05, and 
for T s = .6 and 1.2 for m q = 0.07. 

For the two larger values of T s which we have consid- 
ered, T s — 1.6 and T s — 2.0 the fit of propagators with a 
larger T s in the interval [2:5] are of a poor quality, while 
the fits in the interval [1:5] offer a good description of the 
data, in agreement with results at smaller T s , for t > 2. 
In most cases the propagators for Euclidean time 5 and 6 
where zero within errors, as expected given the large val- 
ues of the mass, and including them does not change the 
quality of the fits. 



m, = 0.05 




T a = .6 


T s = 1.2 


T s = 1.6 


T s = 2.0 


m 


0.85(12) 


0.98(12) 


1.02(12) 


1.05(12) 


x 2 


0.08 


0.08 


0.08 


0.08 


m q = 0.07 




T a = .6 


T s = 1.2 


T s = 1.6 


T s = 2.0 


m 


1.02(8) 


1.7(4) 


2.05(8) 


2.08(8) 


x 2 


1.93 


0.2 


0.2 


0.2 



Table 3. Results for the ++ Glueball at fi — from fits with 
an open vacuum contribution. 



Time interval [2:4] 


m 


T s = .6 
0.93 


T a = 1.2 
0.82 


T s = 1.6 
0.84 


T a = 2.0 
0.89 


Time interval [2:5] 


m 


T s = .6 
0.4(4) 


T s = 1.2 
0.62(25) 


T a = 1.6 
0.68(12) 


T s = 2.0 
0.73(11) 



We summarize the results from the fits at [i = in 
Table[5J For m q = 0.05 , and all T s , we quote results and 
errors from the fits in the [2:4] interval. 

For m q = 0.07, a direct comparison of the results show 
that the glueball masses are indeed heavier than those at 
m q = 0.05. This implies that the propagators at m q = 
0.07 can be compatible with zero already at a distance 4, 
and when this happens is very difficult to obtain reliable 
estimates. For T s — 0.6, 1.2 we quote the results obtained 
by fitting data in the [2,5] interval. They are consistent 
with those in the [3,5] interval, which have however a lower 
significance. For two higher T s values we quote the results 
of the fits in the [1,5] interval (the results in the [2,5] 
interval are m = 2.03(50) , and m — 2.24(65), compatible 
in the large errors), having in mind that these might well 
represent upper bounds rather then asymptotic estimates. 

Finally, we consider also the plaquette correlators eq. 
[T^] where the vacuum contribution was left as a fitting pa- 
rameter. To compare the two fitting forms we considered 
first the same fitting interval [2:4], as well as an extended 
to interval [2:5] , and we show the results for m q — 0.05 
in Table [3l while Fig. [10] illustrates the generic features of 
the results. All in all, subtracting the background allows 
to constrain one parameter fit to zero, thus achieving a 
more precise fit, but the results within the larger errors 
are in substantial agreement. In particular, the fitted vac- 
uum coincides with its direct calculation. 



5.3 Glueballs in the superfluid phase 

We have first performed an effective mass analysis. We 
show in Fig. [IT] a subset of the results, for quark mass 
m = 0.05, and \i = 0.6, and selected values of T s . [i = 0.06 
is well below the saturation threshold, the baryonic den- 
sity is less than one (see e.g. [24]) and the results for the 
amplitudes are far from the quenched ones, which are ap- 
proached only when \x approaches 0.9, as discussed above. 
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Fig. 10. Fit of the plaquette correlator, m q — 0.05, n = 0.0, 
T s = 2.0. The mass value is consistent with the estimate from 
the effective mass analysis and the fits with subtracted back- 
ground and the fitted vacuum coincide with errors with the 
fitted measurement. 



At a variance with the behaviour observed at fi = 0.0, 
here a plateaux begins to emerge, indicating that by in- 
creasing the time distance, and/or the number of smearing 
steps we can approach the asimptotic behaviour represent- 
ing the true value of the mass. 

The fits proven to be rather stable, and indeed we ob- 
served convergence over a rather large interval of T s , the 
convergence to the final results being faster in the time 
interval [1 : 5]. We underscore that for t — 5 the signal 
is very weak, so the interval [1 : 5] is the smallest which 
can be fitted in a meaningful way. We summarise the re- 
sults for the masses over a large set of T s in Fig. [TSJ which 
demonstrates the nice convergence of the results. The ef- 
fective mass analysis give also compatible results, within 
the larger errors. 

Finally, Table [4] collects the results for the glueball 
masses obtained for the fits with /i ^ 0. 

Also in this case we performed fits with an open vac- 
uum contribution, choosing the same interval [1:5] where 
the results of the standard fits proven to be stable with 
increasing smearing parameter T s , and in good agreement 
with effective masses estimates. 

The quality of the fits was very good, however the mass 
results turned out to be larger that those obtained by fit- 
ting the propagator with the subtracted vacuum contribu- 
tion. In addition to this, the fitted vacuum was definitively 
much larger than the direct measurement. Since the fit- 
ted mass was rather large, we should expect not only the 
agreement between fitted vacuum and direct estimate, but 
also an asymptotic behaviour of the propagator reaching a 
plateaux given by such a common value, see Fig. [13l top. 

This is not what we are observing: the plaquette cor- 
relator levels off at a value which is much higher than 
the direct estimate of the vacuum contribution. This be- 
haviour can be explained by a contribution from a light 
mass which is obscured by the background. 



m 




1 2 3 4 5 6 

Euclidean time 



Fig. 11. Glueball effective masses as a function of Euclidean 
time, for /i = 0.6 and m, = 0.05. Results for T a — N 3 x w = 
0,0.1,0.2,0.3,1.2 are plotted slightly displaced one from the 
other, for increasing values of T 3 . 
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Fig. 12. Fit results as a function of the smearing parameter, 
for the two interval considered. The results confirm the ones 
seen in the effective masses, with smaller errors 



To shed light on these observations we performed fur- 
ther fits this time using a constrained value of the mass 
equal to the one coming from the fits with the subracted 
vacuum contribution. In those cases we obtained accettable 
fits, with the correct value of the vacuum contribution, see 
Fig [El bottom. 

A tentative conclusion is thus that it is difficult do 
disentangle the contribution from a low mass excitation 
from an extra contribution to the background. To exhibit 
the correct behaviour we have to constrain the background 
to the measured value. If we do so, either the effective mass 
analysis and the standard fits produce results consistent 
with each other, over a large set of smearing parameters 
and two different time intervals. 

Clearly to confirm this picture we should repeat our 
measurements on lattices more elongated, and with a larger 
size. If the picture outlined here is correct, on a larger lat- 
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Table 4. 

fits 



Glueball masses at fi 7^ extracted from different 
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T s = .6 
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A* = 
A* = 
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II II II 
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Reduced x^ 


associated to the fits above 
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1.73 
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0.82 


0.78 



Table 5. Glueball masses at /i = 0.6 extracted from different 
fits with an open vacuum contribution 
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Ts = 1.2 


Ts = 1.6 


Ts = 2.0 
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1.73(1) 


1.00(5) 


0.89(5) 


0.84(6) 
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1.55(10) 


1.19(7) 


0.89(6) 


0.85(6) 
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tice we should be able to measure a two state signal, and to 
observe with confidence an asymptotic behaviour, stable 
over several smearing steps. 



5.4 The critical region 

As mentioned at the beginning of this Section, the critical 
propagators show a puzzling feature, namely the square of 
the plaquette turns out to be larger than the minimum of 
the propagator itself: if we subract it, as done in the pure 
phases, the propagator would become negative |31j . 

It is rather natural to intepret this behaviour as the 
manifestation of massless modes, which develop at the 
phase transition. In the presence of massless excitations 
there is no reason to expect a plateau in the propagators, 
so the negative region which would emerge subtracting the 
square of the plaquette has no particular significance. 

On the other hand, it is clear that a finite lattice can- 
not accomodate such massless excitations. We only know 
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Fig. 13. The plaquette correlator at m q = 0.05, /1 = 0.6, 
T a — 2.0 with superimposed the fit with an open background. 
The fitted value of the vacuum contribution and its direct mea- 
surement are shown as well, demonstrating a mismatch and an 
overall inconsistence of this fit (top). The same, but with the 
mass value constrained to that measured in the fit with a sub- 
tracted background. The quality of the fit is poorer, but this 
time there are no inconsistencies (bottom) 



that the propagator should become periodical, and it is 
reasonable to assume that it should be continuous with 
all its derivatives at the origin, since there is no discrete 
spectrum dominating the short distance behaviour. These 
considerations suggest to try fits with three parameters 
a, 6, c of the form: 



"X" 



"X" 



"X" 



"X" 



"X 



C(t) =acos(2irt/N t )cosh[b(t - N t /2)] 



(14) 



We show in Fig. Q4] these fits for the raw data, and for 
the smeared ones. The stability of the results is excellent, 
note in particular that the fit in the restricted interval 
[3 : 6] describes perfectly the data at smaller time dis- 
tance, including zero for the smeared results, where the 
ultraviolet fluctuations are suppressed. 

In this picture the glueball propagator in the critical 
region has two complex conjugates poles, whose imaginary 
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Fig. 14. The propagator at m q — 0.05, and /i = 0.2 ~ fi c , with 
superimposed the fits described in the text. The raw (smeared) 
data are shown in the top (bottom) diagram. We obtain excel- 
lent fits to a damped oscillation with period Nt , note in partic- 
ular the smeared results where a fit constrained to the interval 
[3:6] successfully describes the data in the entire interval 



part 2n/Nt goes to zero in the thermodynamic limit. It is 
clear that only a combined finite size and finite Nt scaling 
analysis both in the spacial and in the temporal direction 
might confirm or disprove this intepretation. 



6 Summary and discussions 

We have studied the gluodynamics of two color QCD at 
very low temperature, and nonzero chemical potential. 

We have noted that the Polyakov loop is nearly insen- 
sitive to the transition from the hadronic to the superfluid 
phase, in agreement with the observations of [27], and sug- 
gesting that the superfluid phase shares the same confining 
properties as the hadronic phase in a rather large interval 
of chemical potential. We have not addressed the possi- 
bility of a deconfining transition taking place at an even 
larger fi value and the possibility of a BEC/BCS crossover 
discussed in [271. 



We have observed a gluonic transition - signaled by a 
peak of the plaquette susceptibility - coincident with the 
superfluid transition at fj, c = Beyond this transi- 

tion, in the superfluid phase, the susceptibility reaches a 
lower value than in the normal phase. 

We have studied the scalar glueball propagator in the 
normal and in the superfluid phase. In both phases the 
propagators are amenable to standard hyperbolic cosine 
fits, and related effective mass analysis. In the normal 
phase the results from two different fits - one with a sub- 
tracted vacuum contribution, and the other with the vac- 
uum contribution left as a free parameter - give results in 
good agreement one with another. In the superfluid phase 
the subtraction of the background seems mandatory to ob- 
tain coherent results, and in this case it turns out that the 
lightest excitation in the scalar glueball channel is lighter 
by about a factor two in the superfluid phase than in the 
normal phase. The fit with a free background would pro- 
duce a larger mass, but also an anomalously large fitted 
value for the vacuum contribution, not compatible with 
the direct calculations. 

At the present stage, we can conclude that our results 
indicate a non trivial modification of the gluonic vacuum 
in the low temperature, superfluid phase. 

To put the above considerations on firmer quantitative 
grounds we would need simulations on a larger and more 
elongated lattice, as well as a careful considerations of 
the mixing effects with fermionic channels with the same 
quantum numbers. 

One interesting open question concerns the mass de- 
pendence of the results. We already know that in the large 
chemical potential limit, beyond the saturation threshold, 
we recover the quenched behaviour [25j . The quenched 
limit is also obtained at infinite mass, and in that case 
the gluon fields should not be sensitive to the chemical 
potential. Then, when increasing the mass, the critical 
point would move forward and forward, eventually col- 
lapsing with the saturation threshold at infinite mass, in 
which case all gluonic observables will remain constant 
and equal their fi = value. This suggests the possible 
existence of an endpoint of the transition in the chemical 
potential mass plane, which will be interesting to explore. 

One final comment concerns the interrelation, or lack 
thereof , of these studies with real QCD. Following a large 
7Y C analysis it was proposed [38] that the matter immedi- 
ately above \i c is "quarkyonic" - the authors choses this 
name to indicate a phase with an unusual realization of 
chiral symmetry, which still remains confining. This phase, 
according to the authors, should persist till relatively large 
temperatures, and should not break color symmetry, at 
variance with the ordinary superconducting or CFL phases 
which are indeed features of the very low temperatures, 
and which break color symmetry. The two scenarios of the 
large iV c a quarkonic phase and the small N c superfluid 
phase share a few common aspects: they are confining, 
with unusal realization of chiral symmetry, do not break 
color symmetry and should persist till relatively large tem- 
perature. The physical phase diagram of real QCD at low 
temperatures should in some sense interpolate between 
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the two. It is tempting to regard the superfluid phase as 
a limiting case of the quarkyonic phase of large N c QCD. 
Should this be true, two color QCD might be a better 
guidance to the understaning of the QCD phase diagram 
in the physical case than previously thought. 
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